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SUM  MAR  f 


Solutions  of  the  boundary  layer  equation  for  an  unsteady  flow  have 
previously  ber.i  aotatned  for  only  a  few  boundary  conditions  such  as  those 
v.rucn  exist  in  suddenly  accelerated  or  uniformly  accelerating  flows.  In 
this  paper  a  general  solution  using  the  method  of  successive  approxima¬ 
tions  for  an  arbitrarily  accelerating  flow  is  presented.  The  solution, 
which  is  expressed  in  an  integral  form  including  the  acceleration  as  a 
Chosen  function  of  time,  is  valid  for  both  two-dimensional  and  axially 
symmetrical  flows. 

In  order  to  show  the  feasibility  of  the  solution  one  example  is  pre¬ 
sented  where 

U  (ajt  -  a^t2  *  a}t3  +  a4t4)  h(x) 

i.  e.  ,  the  variation  of  velocity  outside  of  the  boundary  layer  is  a  fourth 
degree  polynomial  in  time  multiplied  by  a  function  h(x). 


V 


LIST  OF  SYMBOLS 

a  =  acceleration  outside  boundary  layer 

h(x)  =  function  depending  on  shape  of  object 

L  -  linear  operator  independent  of  t 

r(x)  -  function  depending  on  shape  of  object 

t  =  time 

1!  =  velocity  outside  boundary  layer 

ut  =  Fourier  sine  transform  ot  u, 

u  =  velocity  component  of  fluid  in  x  direction 

v  -  velocity  component  of  fluid  in  y  direction 

■p*  =  boundary  layer  thickness  as  defined 

r  =  shear  stress 

p  =  dynamic  viscosity 

v  =  kinematic  viscosity 
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INTRODUCTION 


The  boundary  layer  and  continuity  equations  for  a  two-dimensional  flow  in 
the  absence  of  sharp  corners  are 


8u  8u  3u  SU  rT  8U  .  8*u 

aF+ua^  +  vW  =  IT  +ua7+l,e7 


(i) 


and 


8u  8v 
8x  "*  8y 


0 


(2) 


in  which  u  and  v  are  the  velocity  components  of  the  fluid  in  the  x  and  y  direc¬ 
tions  and  v  is  the  kinematic  viscosity.  These  equations  are  to  be  solved  for 
appropriate  boundary  conditions.  Heretofore,  no  exact  solution  has  been  ob¬ 
tained  because  of  the  extreme  difficulty  of  solv.ng  the  nonlinear  boundary  layer 
equation.  Generally,  a  researcher  finds  that  he  must  be  satisfied  witn  a  solu¬ 
tion  involving  successive  approximations.  A  conventional  method  is  outlined  as 
follows : 

(a)  Assume  that  the  velocities  are  the  sum  of  seme  approximation  terms 
u  -  u0  +  u,  +  u*  +  . .  .  .  ,  V  =  v0  +  Vj  +  v2  ►  . .  .  .  (3) 


(b)  Determine  the  values  of  these  approximation  terms  consecutively  from 
the  following  differential  equations. 


8u; 

Ft 

9*u, 

fi(x'V'i 

t). 

L 

-  0,1.2 _ 

(4) 

bx 

8vj 

+  If  = 

0,  t 

=  o, 

1,2.  .  , 

(5) 

fc 

8U 

“  IT 

(6.) 

• 

--fe- 

-  vo 

. 

By 

(6b) 

f* 

da. 

:‘u°  K1 

•a- 

'  u‘eT 

-  vo 

8u, 

* 

8uc 

vl  T3*"  • 

1 8y 

(6c) 

Further  expression*  for  the  functions  fj  can  be  developed  when  higher  order 
approximation  is  required.  The  boundary  conditions  are 

Uo  (x,  0,  t)  =  0,  Ufl  (x,«,  t)  -  U(x,  t),  u,,  (x,  y,  0)  =  0  (7a) 


and 
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Ui(x,  0,  t)  -  0,  ui(x,*,t)  --  0,  ui(x,  y.  0)  -  0  (7b) 


for  i  =  1,2,3.... 

If  the  nonlinear  (convective  acceleration)  terms  in  eq  1  are  small  compared 
to  the  other  terms,  satisfactory  answers  may  be  obtained  by  calculating  the  first 
two  or  three  approximation  terms.  This  may  be  the  case  if  there  is  a  rapidly 

accelerating  flow  where  v  ^73  and  ^  are  large.  Otherwise,  the  calculation  of 

higher  order  approximation  terms  may  be  necessary  to  obtain  the  desired  ac¬ 
curacy 


Solutions  of  the  approximation  terms  have  been  obtained  only  for  a  very 
limited  number  of  problems.  One  of  them  is  the  case  of  a  sudden  acceleration 
of  a  body  starting  at  rest  and  mo  'ing  at  a  constant  velocity  thereafter.  In  this 
case  f0  -  0  for  t  <  0  and  f0  '  U(x)  for  t  >  0.  The  solution  for  u0  is  given  in 
many  of  the  standard  textbooks  on  boundary  layer  theory  (e.  g.  ,  Schlichting, 

I960).  The  solutions  for  u,  and  u2  were  obtained  by  Blasius  (19081  and  Goldstein 
and  Rosenhead  (1936),  respectively. 

Another  problem  is  the  case  of  a  uniform  acceleration  of  a  body  which  was 
initially  at  rest.  The  solutions  for  u0  ,  Uj  and  u^  for  this  case  were  obtained  by 
Blasius  (1908).  More  complicated  cases  involving  some  nonlinear  accelerations, 
where  U  =  h(x)tn  with  n  =  0,  1,2,  3,4,  we  re  studied  by  Goertler  (1944).  Recently 
Watson  (1  955)  extended  the  work  of  Blasius  to  the  cases  where  U(t)  At°  and 
U(t)  --  Aect . 

In  all  the  problems  mentioned  above  the  flow  around  the  body  was  two-dim¬ 
ensional.  Boltze  (1908)  applied  the  method  of  successive  approximations  to  solve 
problems  where  the  flow  is  axiaily  symmetrical.  In  this  case  the  boundary  layer 
equation  remains  in  the  same  fcrm  but  x  and  y  are  the  curvilinear  coordinates 
parallel  and  normal  respectively  to  the  surface  ot  the  body.  The  equation  of 
continuity  changes  to 


in  which  r(x)  is  the  variable  radius  which  specifies  the  con  to  ir  >f  the  body  of 
revolution.  T.ius,  the  continuity  equation  for  the  approximation  terms  is 


which  should  be  considered  instead  of  eq  5  in  solving  the  problem  of  an  axially 
symmetrical  boundary  'ayer.  To  the  writer'9  knowledge,  solutions  for  the  first 
three  approximation  terms  in  an  axially  symmetrical  flow  have  been  obtained 
only  for  the  case  of  a  sudden  acceleration  by  Bolt/.e  (1^08). 

The  writer  has  solved  eq  4  for  an  arbitrary  f  (unction  and  for  an  arbitrary- 
boundary  condition,  namely  an  arbitrary  U(x,  t).  The  only  limitation  is  hat  the 
fluid  should  be  initially  at  rest.  I'he  solution  is  valid  for  both  two-dimensional 
and  axially  symmetrical  boundary  layers. 
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SOLUTION  OF  THE  APPROXIMATION  TERMS 

Before  presenting  the  solution  of  the  first  approximation  term  it  is  neces¬ 
sary  to  give  some  background  information.  Consider  a  linear  differential  equa¬ 
tion  of  the  type 


i  8u 

Lu  r  9f 


(10) 


in  which  L  is  a  linear  operator  independent  of  t.  Let  g(x,  t)  be  a  solution  of  this 

t 

differential  equation.  If  g(x.  0)  =  0,  then  j  F(t-t')  g  (x,  t')dt*  is  also  a  solution, 

o 

This  may  be  proved  as  follows: 

t 

Let  t - 1 ’  =  t".  The  integral  solution  becomes  /  F(t"  )  g  (x,  t-t"  )dt"  .  If  the 

o 

derivative  with  respect  to  t  is  taken, 
t 


A  J  F(l")  g  (x,  t-t")dt"  =  J  F(t")  dt"  .  F (t)g(x. 


0) 


can  be  obtained.  The  last  term  is  zero  since  g(x,0)  -  0.  Since  g(x,  t)  is  a  solu¬ 
tion  of  eq  10 


/  F(t")  dt”  r  /  F(r)Lg(x,  t-t")df ' 


t 

-  L  f  F(t")  g(x,  t-t")dt". 

o 


The  stated  result  follows  immediately.  These  integral  solutions  were  used  by 
Basset  (1888)  to  solve  the  problem  ot  an  accelerating  sphere  in  a  fluid.  They 
are  also  widely  used  in  heat  conduction  problems.  The  function  F(t-t')  is  de¬ 
termined  from  the  boundary  conditions. 

Now  the  problem  cf  the  unsteady  fluid  flow  caused  by  an  arbitrary  and 
rectilinear  motion  of  an  infir  itely  long  flat  plate  will  be  investigated.  The 
Navier-Stokes  equations  reouce  to 


lu  8*u 


(ii) 


in  which  u  is  the  velocity  paralle’  to  the  plate.  If  the  plate  is  accelerated  sud¬ 
denly  from  rest  and  its  velocity  is  kept  at  unity  thereafter,  the  boundary  con¬ 
ditions  are  u(y,  0)  -  0,  u(0,  t)  -  1  and  u («*  t)  -  0.  The  solution  which  can  be 
found  in  many  standard  textbooks  (e.g,  Schlichting,  196£)  is 


u  1  -  erf  e| 


(U) 


I 


'X 
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where  rj  =  £y(vt)'2  .  For  an  arbitrary  acceleration  the  boundary  conditions  are 
u(y,  0)  =  0,  u(0,  t)  =  U(t)  and  u(oo,  t)  =  0.  The  integral  solution  is 

t 

u(y,  t)  =  J F(t-t')  [  1  -erfq’jdt' 
o 

where  r|'  =  jy( vt')_ a  ,  According  to  the  second  boundary  condition, 
u{  0,  t)  -  /  F(t-t')dt'  =  U(t). 


Thus, 

F(0  =  =  4t) 

where  a  is  the  acceleration.  The  solution  becomes 


u(y.  t) 


j  (1  -  erf  r) 1 )  a (h— t*)dt* 
o 


where 


erf  t|* 


2lT~2 


(13) 


Th;s  result  will  be  used  to  obtain  the  solution  of  the  first  approximation 
term.  The  differential  equation  for  the  first  approximation  term  is 


9u0  82u0  8U 

it  '  =  TT 


and  the  boundary  conditions  are  u0  (x,  0,  t)  =  0,  u0  (x,  w,  t)  =  U(x,  t)  and 

u0  (x.  y,  0)  =  0.  By  introducing  a  new  variable,  u0'(x,  y,  t)  =  u0  (x,  y,  t)  -  U(x,  t), 

thie  differential  equation  reduces  to 


(15) 


and  the  boundary  conditions  change  to  ^'{x,  0,  t)  =  -U(x,  t),  u0'(x,  °o,  t)  =  0  and 

Uo'(x,  y.  0)  =  0.  Note  that  U(x,  0)  =  0. 

Equations  11  and  15  and  their  boundary  conditions  are  the  same  except  that 
the  second  boundary  condition  of  cq  15  has  a  minus  sign.  Thus  the  solution  for 
the  first  approximation  term  can  be  readily  obtained  in  the  form 


u<> 


J  -ftTi  erf  |y[v(t-t")] dt". 


(16) 
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Inorder  to  calculate  higher  order  terms  of  eq  I  a  Fourier  sine  trans¬ 
form  with  respect  to  y  is  used.  The  Fourier  sine  transform  and  its  inverse 
transform  are  respectively 


Z  10 

J  Uj  sinysdy  and  u  [z~  J  In  sinysds. 


(1?) 


Multiplying  eo  4  by  sinys,  integrating  the  second  term  twice  and  using  the 
boundary  conditions  on  uj(x,y,t)  gives 

dC  °° 

-y1  +  vszu{  =  /  fj(x,y,  t)  sinysdy.  (18) 

o 

Inis  is  a  first  order  non-homogeneous  Linear  differential  equation,  with  solu¬ 
tion 


t  ao 

TT  i  =  exp(-vs2t)  J  exp{vs2t')  J  £■  (x,y,  t')sin  ys  dy  dt*  +C  exp(-vs2t). 
o  o 

According  to  the  third  boundary  condition,  when  t  =  0,  uj  =  0  and  therefore, 

Hi  =  0.  This  gives  C  =  0.  Using  the  inverse  Fourier  sine  transform  the  solution 
u-  is  obtained: 

uj  =  ~  Jj  J  exp[  -  v  s2(t-t’)]  fi  (x,  y',  t')sin  ys  sin  y's  dy' dt' ds.  (20) 
ooo 


Integrating  with  respect  to  s  this  equation  reduces  to 


SKIN  FRICTION 

The  skin  friction  can  also  be  expressed  in  series  form: 


(ZZ) 


in  which  r  is  the  shear  stress  and  p  denotes  dynamic  viscosity.  The  derivatives 
of  the  approximation  terms  at  y  =  0  are 


1 _ 


I 


0U 

(t-t"/2 


dt" 


(23) 
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and  for  i  =  1,2,3... 


t 


Substituting  expressions  23  and  24  for  uj  in  eq  22, 


(24) 


dy' d  t* 


(25) 


is  obtained.  This  is  a  series  solution  for  the  skin  friction.  The  number  of  ap¬ 
proximation  terms  depends  or.  the  type  of  motion,  the  shape  of  the  body  and  the 
degree  of  accuracy  required  by  the  problem. 


APPLICATION 


U  =  (a|t  +  a2t2  +  a^t3  +  aAt4)  h(x) 

In  this  example  the  terms  for  i  =  0,  1  and  2  of  series  7  are  calculated.  The 
first  integral  term  is 


t 

I 


3U  +  u  au 

— — -5^  dt'  =  {  2a,  +  2.  667  a2t  +  1 . 067  (3a,  +  --  a,*)t* 
(t-t  ')*  '  dx 


0.914  (4a*  +  2^|  aja2)tJ  +  0.813  (2a,a,  +  a**)  ^  t4  +  1.478(aza,  + 


+  aia4^t’  +  °-682  (Za2a4  +  as*)^t4  +  1>2?3  a^J-t7  + 

+  0.599  .  (26) 

The  zeroth  order  approximation  term  for  velocity  is 


t 

u0  =  h  f  erf - ^ - r  (aj  +  2azt"  +  3a,t"*  +  4a4t"*)  dt".  (27) 

Z[v  (t - 1‘ 1 )] a 

0 
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The  integration  leads  tc 

u0  =  h{2ta,  F0(r))  +  4a2t2F,(T>)  +  6aJt5F2(T1)  +  Sa^Fjfa)} 

where  F0  ,  Fj,  F2  and  Fj  are  functions  of  q  -  — ^ — r-  , 

2(vt)2 

F o  =  erf  Tj  (r|2  +  j)  +  tr*2  t}  e  11  -if 
Fj  =  erf  q  (1-  +  if  +-)+  w'a  +  yrf)  -  rf  -  i-tf 

F2  =  erfi^i+rf +  |rf +  ~rf) +w  2eV  (j-|T1+|!rj»  + 


45 


4  5\  2  ^4  4  i 

45  ’l  )  *  11  '  3' n  ’  45^ 


Fj  -  erfr|  (^  +  rf  +  V  +  Js^  +  +  "  2  (rifi  '1  + 


37 
42 1 


140 

l  8 


(28) 


(29a) 

(29b) 


(29c) 


(29d) 


ri(  2 

The  next  step  is  the  calculation  of  3-f .  To  do  this  expressions  for  ■?-%-  , 
a  ,  .  9y  dYdx 

Vq  and  Q  ;u  are  needed. 

It  can  be  shown  that 


92u0 

9y9i 

-  r  (v  0"  *-  ^  {a1G0t+ 2a2G,t2  +  3a3G2tJ4  4a4Gjt4> 

(30) 

in  which  G0  , 

Gj,  G2  and  Gj  are  some  functions  of  t). 

G0 

.1  > 

=  erfr|  +  2tt  2  e’1"1  -  2r| 

(31a) 

G, 

=  erf  r)(yr)*  +  2rj)  ■>  it  2  e'^2  (|+|if)  -  |rf  -  Zq 

(31b) 

g2 

r  erftl(-^T19+  |rf  +  2^)  +  it  2  e"H2  (||  +  ~r)2  4-^r!4)  - 

?  ®  j  8  5 

2t>  '  3 ' n  '  15^ 

(31c) 

G, 

=  erfTl(li^Tl7  +f|rl5  +  4t1S  +  2rl)  +  "  2  e"'!2  (jf  +  '-jT'S2  + 

32  4  16  6\  5  a  1  72  §  16  7 

+  IT’1  +I05^  ■  2ti  •  4t1  *  45^  •  Tol^ 

(31d) 

and 


•& 


| 


I 

V 

f- 

■f 


.1, 

I 

.  +• 


—  »•>  -Mb 
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02  ^  ^ 

'  TZ  ^a‘s°  *  2a*  ^s°  ‘  si)f  *  3aJ  (so  'Zsi  *  s2)t2  * 


4a4(s0  -  3s j  *  3s2  -  s,)  }.  (32) 

In  which  s0  ,  s,,  s2  and  s3  are  some  functions  of  q. 

s0  =  -2erfcq  (33a) 

.I  2 

st  =  W  erf c  q  -  4 IT  2q  e'9  (33b) 

s2  =  -|q4erfcq  +  (2q*  -  (33c) 

s3  =  erfcri  -  itt**  (fri5  -  §nJ  +  n)e  *12.  (33d) 


Using  the  continuity  equation  (eq  5)  v0  can  be  calculated. 

vo  =  -4(v  t)1  ~  (aiV0  t  +  2a2(V0  -  V,)t2  +  3a3(V0  -  2  V,  +  V2)t3 

+  4a4(V0  -  3  V,  +  3V2  -V,)t*  }  (34) 

in  which  V0  ,  V,,  V2  and  V3  are  some  functions  of  r). 

V0  =  erfq(iq+iq3)  +  rr'*  (1  +  ±q2)  e‘^  -  jq3  -  }  (35a) 

V,  =  erfqfiq-^qWM^q2  -i^e^ 


1  s  1-4- 
+  Ts1'  '  571 


(35b) 


4 

315 


7  .  1  -* 

•q  tT  7t  ' 


(35c) 


Vj  =  erfq  (fn  ‘955Tl’)+Tr"2  (d  +  TF7 


>  ,l  ■  1  2  -  * - 4  ,  1 


9  +  252^  *  ^30 13  ’945^ 


^  8  v  - -p^  Z  a  1  -4 

-  9459  )«  71  +  945^’  -  9  " 


945 


(35d) 


di 


N°w  ^  can  be  calculated: 


'  ~8y  =  +  4^3,32X2  t-  6t3a,a3X3  +  8t3a22\4  +  8t4a,a4\5  + 

1 

+  12f4a2a3\6  +  1  6t*a2a4X7+  18t5a32\8  +  24t‘a3a4X,  +  32t7a42^j0}h^  (+)* 

(36) 
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in  which  \j,  i-1  to  10,  are  some  functions  of  q  which  involve  Fj,  Sj,  Gj  and  Vj, 
j  =  0  to  3. 

2 

3f  r  v'  1  -i 

The  next  step  is  to  multiply  with  exp  j-  4y(t.ttyj  2  anc^  to  2nte* 

grate  with  respect  to  t'.  Substituting  t/(t-t')  -  a  and  carrying  out  the  integra¬ 
tion  the  following  result  is  obtained. 


t  00 


o  o 


exp 


[•*£d 


(t-f)2 


3fj 

3y 


dy'dt'  =  -|4t2aj2I,  +  8t*a|a2l2  + 


+  12t4a|a3Ij  +  I6t4a22l4  +  l6t*a1a4Is  +  Z4tsa2aJI6  +  321^3*1 7  + 

+  36t6as2I8  +  48t7a,a*I,  +  64tl42Il0  |  h  t2  (37a) 


where 


oo 

*  ■  I  *. 


<j>  dr|  for  i  =  1 , 2,  ....  1 0 


(37b) 


where  i  refers  to  the  terms  of  the  series  in  eq  37a  and  <j>  is  a  function  of  q  cal¬ 
culated  from  the  integral 


4>(j.ri)  =  e1"!2 


(38) 


where  j  is  the  numerator  of  the  power  of  t  in  eq  36,  every  term  of  which  can  be 

expressed  in  the  form  of  aman\jt^^.  The  functions  $(j,q )  are  tabulated  in 
Appendix  A. 

The  integrals  Ij  asymptotically  approach  constant  values  as  q  —  °ct  It  can 
be  shown  that  if  they  are  integrated  from  zero  to  properly  chosen  values  of  q*, 
the  values  of  the  integrals  remain  practically  unchanged.  For  example 

I,  =  0. 0536  for  q*  —  oo  I,  =  0.  0534  for  q,*  =  1.44  (39a) 

I2  =  0.0429  for  q*  — »  I2  =  0.  0428  for  q,*  =  1.44  (39b) 

I*  =  0.  00952  for  q*  —oo  I,  =  0.  00951  for q2*  =  1.23.  (39c) 

The  values  of  q2*  and  q2*  refer  to  the  boundary  layer  thicknesses  for  the 
cases  where  Uj  =  hatt  and  U2  =  ha2t2  respectively.  They  are  calculated  from  the 
relationships  u01/U,  =  0,99  and  uM/U2  -  0.99  ana  used  as  upper  limits  for  the 
first  and  fourth  integrals  since  X.2  <f>,  and  involve  only  the  functions  due  to 

the  terms  linear  and  quadratic  in  t  respectively.  The  second  integral  involves 
the  cross  product  terms  and  therefore  the  higher  value  of  q*  is  chosen. 


1 

?• 
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The  values  of  the  other  integrals  are 


0.0237  for  qt*  = 

1.44 

(39d) 

0.  0151  for  q,*  = 

1-44 

v39e) 

0.0112  for  q2*  = 

1.23 

(39  f) 

0.  00751  for  q2*  = 

1.23 

(39g) 

0.  00345  for  q2*  = 

1. 09 

(39h) 

0.  00477  for  q,*  = 

1. 09 

(39  i) 

0.  00168  for  q4*  = 

0.  99. 

(39  j) 

Now  the  next  term  of  eq  25  will  be  calculated.  The  first  order  approxima¬ 
tion  term  for  velocity  should  be  calculated  by  using  eq  21.  It  can  be  shown  that 


=  +  2aja2S2  (q)tJ  +  2a1a,S,  (ri)t4 


+  a^S^Jt4  +  2a1a4S5(q)t*  +  2a2a3S4  (q)t*  +  2a2a4S7  (q)tfc 

+  a3ZSd  (r|)t^  +  2aja4S9  (q)t7  +  a^Sjj  (q)t  }  (40) 

where  Si(q),  i  =  1  to  10,  are  some  functions  of  q  calculated  from  eq  6b,  The 
terms  of  fj  are  to  be  multiplied  with  the  kernel  in  eq  21  and  integrated.  Since 
the  Sj(q)  functions  are  quite  involved,  they  are  approximated  to  fifth  degree 
polynomials  in  q  by  using  the  least  squares  method. 

The  next  step  is  the  integration  in  accordance  with  eq  21.  Each  term  has 
the  form  of 


t  oo  ,mi  y>  .  n 

UIL  =  hg  ‘  f  - {K}  dy'df  (41) 

oo  '  ' 

in  which  K  denotes  the  kernel.  Changing  the  variables  t’  -  a2t  and  obi  Irving 

i  I 

that  y'  =  2{v  t')2  q1  and  y  -  2(vt)?q,  u  j  j  can  be  written  in  the  following  form 


0  .  . 


(42) 


Since  the  value  of  u2j  becomes  small  at  the  edge  of  the  boundary  layer,  the 
upper  limit  of  the  integral  is  changed  from  oo  to  qj*  and  the  integration  with  re¬ 
spect  to  q'  is  performed  for  each  term  of  the  polynomial,  i.e.,  for  the  values 
of  n  varying  from  0  to  5.  The  results  are  shown  in  Appendix  B.  Next,  the  in¬ 
tegration  with  respect  to  a  is  carried  out  by  a  digital  computer  at  chosen  values 
of  q  with  the  values  of  m;  varying  from  2  to  8.  Thus,  eq  42  reduces  to 
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,  -4-  ,  dh  m:rl 

uJi  =  2  hr-  t  i 


5 

£  A 
0 


in  ^ni  M 


where  the  values  of  Cn^  are  tabulated  in  Appendix  C.  The  values  of  Uj-  are  cal¬ 
culated  at  chosen  values  of  r)  and  the  resulis  are  approximated  to  tenth  degree 
polynomials  in  q  by  using  the  least  squares  method.  Thus,  an  approximate 
expression  for  u,  is  obtained. 

U‘  =  JT  WgiOl)1*  +  Za^gjMt4  +  ^ajajgj^Jt* 

+  +  ^aia4gj  (h)t6  +  2aiaJgi (q)t*  +  2a2a4g7(q)t7  +  a,zgs(q)t7 

+  2a,a4g,(r1)t8  +  a42g10  (q)t’ }  (44) 

where  gi(r) ),  i  =  1  to  10,  are  tenth  degree  polynomials  in  q.  Since  in  this  cal¬ 
culation  process  two  approximations  with  the  least  Squares  method  are  involved, 
it  is  advisable  to  check  the  accuracy  cf  the  expression  for  uj.  To  do  this,  Uj  is 
inserted  in  eq  4  and  it  is  found  that  for  all  practical  purposes  the  equality  is 
verified.  In  this  connection  it  should  also  be  pointed  out  that  the  values  of 
gj(i|)  should  be  the  same  as  the  values  of  Blasius'  function  multiplied  by 

jtt2  .  In  order  to  make  the  comparison  both  of  these  functions  are  plotted  in 

Figure  1.  The  values  of  the  functions  up  to  q  =  0.  8  are  practically  the  same 
and  the  curves  deviate  slightly  from  each  other  above  q  =  0. 8.  This  is  further 
evidence  that  the  polynomial  approximations  are  satisfactory.  Incidentally  the 
value  of  Blasius'  function  at  q  =  1 . 0  is  0.  069  and  not  0.  02  as  tabulated  by  Blasius. 
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~dy  can  calculated 


Ty  “  £  (r.^1  jia.'t*  (r.^Sl  ♦  s0  e,  |  +  4a1ia1t*/1 

|<Fo3^  *■  <FI  ^  ’  sue,)|  •  4a,*a,t,lA  |(F0  ^  s0e,)  ,  1.5 

(F‘  $  *  satei ) J  •  j(F„  i£*  +  soe*)  *  4(F, 

*  WV"'‘  {^0  Iff  -*.4,)  -2  (F,^  .  slte,)|  .  4a,a,a,fV* 

l^F»  d5f  '  *«  e‘>  *  2(FI  ^  -  ‘itej)  »  3(F;  iSi  *  s2tej)j  ,  4a,’t'>/i 

iF‘^  I<f»3^  *»,«,)*  3(F,iSi  ,  ,itej) 

*  4<F>l?f*  [(Ft^.5(e,).6(F2iSa  .  ,!tej)| 

*  l(F'  1%  *»•”  (F,^  .  4.lW„/« 

((F»1^  *s«-l)  *  3(Fil^  *  s„e5)  *  4(F,i8a  ,  S)te,)|  ,  8a,  V"'* 

l(F‘^h  *  *>'*’>  *  (F>^‘  *  •..*.)!  *  4a,a,‘t*’/i  [(F,|S1  ,  ,lieJ 

*  3<F<^  *  ««e.)|  *8.,.,V'V‘  |(F,iS»  l.5(F,|tl  a,lte 

*  «*>!$  *  ».,••»/*  |r,|a  . 

I(F*  *  8<Fy^  •  •  4a,a,*(“/»  j(F,2|Su.  ,  Sii,i() 

*  “,V"*  |(F.^  •  -a,.,.)  •  f(F,^  ,  8a,.,W, 

}i.,V/-|c:g,  -  v„^l} 
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+  1  (Go  g2  -  Vo  yjt  )  +  (G,g, -Vlt^)|+  4a,*a,t»/* 

(Go  gj  *  V0  )  +  1 .  5(G2gl  -  Vn^  )]  4  2a1ai2t1,/i  j  (G0  g*  -  V0  ^  ) 
+  4(Glgl-V,t^)l  +  4a,2a4t»A  f(G0gs-vJ&)4  Z(Gjgl  -  V,t  ^ )1 


4aja2a5t,J/2  (G0  g4  -  V0-r-j^)  +  ^(Glgj  -  Vlt  -rfy )  +  3(Gzg2  -  V2t^-^-)  ! 


4a1a2a4tis^2  (G0  gl  -  V0-T-$r  )  +  Z(G,gs  -  +  4(G,g2  -  V,t^r%)l 


+  Zaia,2tis/2  |^(G0  g8  -  V0  ^ )  4-  6(G2g,  -  V2t  |  4  4a1a,a4t*7/2 

(G0g,- V0^)  +  3(G2g5-V2t^)  4  4(Gjgj  -  VJt^)l  4  Za.a^t"/2 

(G0S10  -V0^L)  4  8(G,gs-VJt^)l  4  4a2W  [clg4-Vlt^  1 
4  8a22a,t15/2  UGlgt-  Vlt^)  4  0.75  (G2g4-V2t^)l  4  8a22a«t>7/2 


[(gi8i‘ Vjt^1)  +  (Gjg4-V,i^)1  4  4a2aJ2t11/1  (G,gg  -  V,t ~t4* ) 


‘t  dq*  '  '  164  *  dqz  '  \ 

+  3(G,g6- V2t^)  !  4  8a2aJa4t1,/z  J^(G,g4  -  V,,^)  4  1 .  5(G2g7  -  V2t  ) 
+  Z(C’J8*'  +  4a2a42t2'/2  |^(G,g10  -  V,t^*)  4  4(G,g7-  V,t^)] 

7  6a3J‘‘,/j  [c,g,-  Va^]  4  1  ZaJza4ti*'2  ^(C2g,-  Vit^  )  +  j 

(GjB.  •  V,,  t-41  )"j  4  '  a.a^t21''2  ^(G,g,?  -  V,.  4  |  (Gjg,  -  V..  }] 

*  |c,gle  -VJt^]J  (45) 


where  e 


i  :  ~  ^  8idn.  *U  :  »0  -  *1. 


—  ZSi  t  8 


i  ♦  81.  »jt  ;  30  •  }»i  ♦  3»*  - 


vit  =  Vo  -  V„  Vlt  --  V0  -  ZV,  4  V2,  v,t  :  Vn  -  3V,  4  3 V2  -  V,. 

Each  term  of  this  equation  has  the  form  of  am  an  a*  t^  Xp(q)  where  X  is 
the  function  of  q  in  the  brackets.  Each  term  is  to  be  multiplied  with  ^ 
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f  VJ  ,  ,  .1 

exp  [  4w(t-t')  ]  ^'ll)  2  and  integrated  with  respect  to  i‘.  By  substituting  4~t  =  a 

and  integrating  with  respect  to  a,  o(j,n)  functions,  where  j  is  the  numerator  of 
the  power  of  t,  are  obtained.  They  are  tabulated  ,n  Appendix  A.  Thus, Tach 

term  can  be  wrill»n  ir.  farm  of  am  a„  ak  0*2  ~f  *(j,n)  Xp  (n)d,.  Changing  , he 

V  “d  •"«  -egrauan,  ,he 


!"FT??rWl  st.  Mr 

— dy  d''  =  -dhti  )*.*•*  [o. 


00956  h- 


+  0.0282  (dxjfJ  +  a»  a^t5  [o.0229  hl£  +  0.  0828  (3^/]  +  a,2a,t6  [o.  01  95  h  ^ 
+  °-°531(^)1  +  ».».*t‘[0.0182hg  ♦  0.0508  0]  ,  a;‘a4t’[O.Ol72h0 
+  °-°455©2]  +  a,W’[O.O3O9h0  ♦  0.08560]  *  a,V 

[°-  00485  +  °‘C1360]+  atW  [°.  0266  h ,  0.  108(g)*]  ♦ 

a,aJtf[O.O134h0+  0.0405  0],  a^f  [o.  0 11  3  h  g  *  0 .  0649  0/]  * 
ala,a,t’[o.02s8hi^  *  0.  0677  0],  a2‘a4t»  [0.  0 1 1  0  h  0 , 0.032  0], 

a*ai*t*[o.  0*  56h^  ,  0.0311/^)'].  a2a,a4t*°  [o.O24nh0 . 0.0662  0], 

a,-V°Jo.00  3s7hl£  >  O.O°n430],  a,a4V»  fo.  0 1  4  h  0 . 0.0374  0*]. 
a2a4it"fo,0.  0  ,  0.04320]  .  a,V"  [0.  0 1  2 .  10  .  0.03030']  . 
Wtu[o.OHthg  .  0.0  3420*].  *,»i»[o.OO57h0  .  O.PIH0](. 


Higher  order  approx ini.it urn  terms  for  ▼  and  u  c.vi  be  ca. dilated  m  ;he 
san.e  n  tier.  Po.vnomia.  app ;  ox  in:  at ,  or;  s  are  quite  accurate  and  can  be  sue  - 
'  d^h  /dX  !  aPPr-”“^ation  terms.  Gene  rail  v  the  coefficients  of 

^  dx*  (cTx)  r'4Pl^>^  approach  .heir  tma.  values.  For  example,  the  first  set 

•!  C  neff.t  lent »  changes  only  bv  between  n  0.  8  and  n  1.4.  Therefore  ,|,ght 
is.  re  pans  le  s  in  u  or  o.igh  val  ies  ot  such  as  the  discrepancy  ab.-.e  n  0.8  in 
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the  expression  for  g,(n),  nave  very  little  effect  in  the  final  answers.  However, 
if  more  accuracy  is  required,  different  ranges  of  q  may  be  considered  and 
separate  approximations  for  these  ranges  may  be  made.  Incidentally  the  coef¬ 
ficients  of  b  and  for  uniform  acceleration  should  be  the  same  as  tnose 

given  by  Blasius.  B!as;us'  values  are  slightly  different;  they  are  0.  010  and 
0.026  respectively.  Apparently  Blasius  has  made  a  slight  numerical  error  be¬ 
cause,  when  these  coefficients  are  calculated  according  to  Blasius1  own  formula, 
the  sam"  values  as  those  presented  in  this  paper,  are  obtained. 


CONCLUSION 

The  integral  solutions  (eq  16,  21  and  25)  are  useful  in  calculating  the  suc¬ 
cessive  approximation  terms  for  the  shear  stress  in  an  unsteady  laminar  boun¬ 
dary  lay  'r.  These  integral  solutions  are  valid  for  arbitrary  accelerations  and 
in  this  paper  the  i' 'eg  rations  are  carried  out  for  the  first  two  successive  ap¬ 
proximation  terms  in  an  example.  The  example  is  chaser  so  that  it  mav  be 
useful  to  the  readers  in  approximation  of  the  other  velocity  variations. 

Generally,  it  is  difficult  to  calcula'  '  ese  integrals  without  any  approxi¬ 
mation.  However,  as  shown  in  this  paper,  polynomial  approximations  leading 
to  accurate  results,  for  most  practical  purposes,  can  be  made.  Higher  crder 
successive  approximation  terms  can  be  readilv  obtained  in  a  similar  manner. 
Since  these  calculations  can  be  carried  out  with  a  high  speed  digital  computer, 
the  method  may  prove  quite  useful  in  calculation  of  the  successive  approxima¬ 
tion  terms  in  unsteady  boundary  layer  problems  where  separation  does  not  take 
place.  Higher  order  successive  approximation  terms  may  no!  oe  needed  for 
the  cases  where  separation  takes  place  since  generally  reparation  occurs  a 
short  time  after  'he  onset  of  the  motion. 
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$  F unctions 

4.(3,  n)  =  Z,  '  2Zj  +  Z7 

-}>{5,  r,)  =  Z,  -  3Z,  +  3Z7  -  Z, 

4>(7,t,)  =  Zs  -  4Zj  +  6Z7  -  4Z,  +  Zj, 

4.(9,  t,  )  =  Z,  "  5 Zs  +  1  0Z7  “  1  OZ,  +  5Z,j  ~  Z,j 

<(»(  1 1 , T^)  =  Z,  -  6Z,  +  1  5Z7  -  20Z,  +  15ZU  -  6Z,,  +  Z15 

4>(  1  3,  t)  -  Z,  -  7Zj  +  2  1  Z7  -35Z,  +  35Z„  -  21Z„  +  7Z1S  -  Z17 

*(15,1,)=  Z,  -  8Z}  +  28Z7  -  56Z,  +  70Zn  -  56Zn  +  28Z1S  -  8Z„  +  Z„ 

4.(17,  r,)  =  Z,  -  9Zj  +  36Z7  -  84Z,  +  126Zn  -  126ZU  +  84Z15  -  36Z17  +  9^„  *  Z2, 

4>(19,n)=  Z,  -  10Z,  +  45Z7  -  120Z,  +  210Z„  -  252Zn  +  210Z1S  -  l20Zn  +  45Z„ 

-  1  cz21  +  Z2J 

4.(21, n)-  Z,  -  11Z,  +  55Z-,  -  165Z,  +  330Zn  -  462Za  +  462Z1S  -  330Z„  +  165Z„ 
-5^Z2|  +  1  1  Z2,  -  Z2S 

4.(23,  t,)  -  Z,  -  1.2Zj  4-  66Z>  -  220Z,  +  495Zn  -  792Z„  +  924Z1S  -  792Z1?  +  495Z„ 
-22 0Z2,  +  66Z2i  -  1 2Z25  +  Z27 

4>(25,t,}  =  Z,  -  1  3ZS  +  78Z7  -  286Z,  +  715Z„  -  1287Z,,  +  1716ZIS  -  1716Z,7 
+  1287Z-,  -  7:3Z21  f  286Z2J  -  78Z25  +  13Zr  *  Z„ 

All  Zj  are  functions  of  r,  as  defined  below: 

Z,  =  2.0  -  3.  54490770r,  e^  erfc7, 

Zs  =  0.  66566667  -  1.  3 333333  Jr,2  +  2.  363271  79r)J  er'?  erfc  r, 

Z?  =  0.40  -  0.2666666  V  +  0.  53333333t,4  -  0.  94530871r|5  s’1  erfc  n 

Z,  =  0.28571428  -  0.  i  142857 V  +  0.  0761 9048rj4  -  0.  i  523809 V  4-  0.  27008826t11 

e*1  erfcn 

Z„  =  0,22222222  -  0.  0634920 V  +  0.  02539682tf  -  Q.  01  693 122*,4  +  0.  03386243t,»- 
0.  06001  erfc  r, 

Zj,  =  0.  18181818  -  0.0404040V  +  0.  01 1  54401n*  -  0.  00461 76V  +  0.  Q0307840r)»  • 
0.  0061568  V0  +  0.0109l266r)ue’1*  erfcr, 

7,s  =  0.15384615  -  0.02797203!,*  +  0.  0062 1  fe01r|4  -  0.  001  77600r,*  +  0.  00071  040n»- 
0.  00047353!,'°  +  0.  0009472  V*  *  0,  001 67383r,‘1  e1^  erfc  r, 

Z,T  =  0.13333333  -  0.  02051282r,2  +  0.0037296V  -  0.  0008288V  +  0.  0002368V  ~ 
0.  00009472r,io  +  0.  000063 1  5t,u 

Z„  =  0.11764706  -  0.01568627^  -  0.  00241  327r,4  -  0.  00(3387  V  +  0.  00009751r,» 
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